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Abstract 

The total canonical (Liouville-Henneaux-Gibbons-Hawking-Stewart) mea- 
sure is finite for completely nonsingular Friedmann-Lemaitre- Robertson- Walker 
classical universes with a minimally coupled massive scalar field and a positive 
cosmological constant. For a cosmological constant very small in units of the 

m 

. square of the scalar field mass, most of the measure is for nearly de Sitter so- 

lutions with no inflation at a much more rapid rate. However, if one restricts 
to solutions in which the scalar field energy density is ever more than twice 
the equivalent energy density of the cosmological constant, then the number 
of e-folds of rapid inflation must be large, and the fraction of the measure is 
low in which the spatial curvature is comparable to the cosmological constant 
at the time when it is comparable to the energy density of the scalar field. 

The measure for such classical FLRWA-</) models with both a big bang and 
a big crunch is also finite. Only the solutions with a big bang that expand 
forever, or the time-reversed ones that contract from infinity to a big crunch, 
have infinite measure. 
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1 Introduction 



Starting with the Louiville measure, the procedure of Henneaux [T] and, in more de- 
tail, of Gibbons, Hawking, and Stewart [2], provides a natural canonical measure on 
the set of classical universes. For a minisuperspace model of a A; = +1 Friedmann- 
Lemaitre-Robertson- Walker (FLRW) geometry (homogeneous and isotropic three- 
sphere spatial sections) with a time-dependent scale size a{t) and a single minimally 
coupled homogeneous massive scalar field the total measure is infinite, and 

Hawking and I showed [3] that all but a finite measure have arbitrarily small spatial 
curvature (or arbitrarily large size a) at any fixed positive value of the energy 
density of the scalar field. However, we also showed [3] that both the set of inflation- 
ary solutions and the set of noninfiationary solutions have infinite measure, with an 
ambiguous ratio: the Louiville-Henneaux-Gibbons-Hawking-Stewart canonical clas- 
sical measure gives an ambiguous prediction for the probability of inflation. 

Gibbons and Turok ^ sought to remove this ambiguity by identifying universes 
in which the spatial curvature is too small to be distinguished, though in 

another way of looking at it, it is not clear why one should be justified in identifying 
universes that have scale sizes a (and hence also spatial volumes that are proportional 
to a^) that are so large and different. Indeed, Turok has abandoned this approach 
and is working on another [5]. In both of these procedures, one gets a classical 
probability of inflation that goes approximately inversely with the volume expansion 
factor during inflation, i.e., roughly as exp (— 3A^), where is the number of e-folds 
of inflation glE]. 

Here a very simple alternative restriction of the set of classical solutions is consid- 
ered that does not depend on any identifications of what might not be observationally 
distinguished and does not depend on any arbitrary choice of finite ranges. In partic- 
ular, I consider the set of classical FLRW solutions that are completely nonsingular, 
with neither a big bang nor a big crunch. For FLRW solutions with a homogeneous 
minimally coupled massive scalar field and with a cosmological constant, nonsin- 
gular solutions with positive canonical measure occur only for the k = +1 FLRW 
geometries (allowing a{t) to have a minimum value) and for a positive cosmologi- 
cal constant A (allowing the universe to expand forever in both directions of time). 
(There is an uncountable set of perpetually bouncing solutions for a homogeneous 
massive scalar field minimally coupled to a = +1 FLRW geometry with A = [6], 
but this apparently fractal set is discrete and so has zero canonical measure.) This 
restriction of positive measure for nonsingular solutions to = +1 and A > is also 
true for any homogeneous minimally coupled scalar field with a canonical kinetic 
term and a potential term with one single extremum that is a minimum of zero 
value, but here we shall focus on the homogeneous free massive inflaton scalar field 
(f){t) of mass m, taking the cosmological constant to be A = 3/6^ = 3m^A > with 
a characteristic length scale b and a dimensionless rescaled cosmological constant A. 

This paper shows that the set of such nonsingular classical FLRW universes has 
finite canonical measure. Therefore, if one restricts to such cosmologies with neither 
a big bang nor a big crunch, one can get unambiguous finite fractions for any subsets 
one chooses. 
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Friedmann-Lemaitre- Robertson- Walker closed model 
with massive scalar field and cosmological con- 
stant 



The Friedmann-Lemaitre-Robertson- Walker (FLRW) closed model with massive 
scalar field and possibly also with a cosmological constant has been analyzed many 
times previously [71[Hll9l|6l[ini[IIl[12l[l3l[Tl[T5l[T6], so many of the dynamical 
equations I shall give below have been previously given, along with much of the 
qualitative behavior. 

The k = +1 FLRW spacetime metric is 



-N\t)de + A\t)dVtl = -±-{-n\t)dt' + a''{t)da 



I), 



(1) 



where N{t) is the lapse function, A{t) is the physical scale size (which is what is 
usually called a{t), but I am reserving that for the rescaled scale size), dfl^ is the 
metric on a unit 3-sphere that has volume 27r^, n{t) = mN{t) is a rescaled lapse 
function that is dimensionless if t is taken to be dimensionless, and a{t) = mA{t) is 
a rescaled scale size that is also dimensionless. 

Using units in which Ti = c = 1, but writing Newton's constant G = rripi or the 



Planck mass mpi = G"^/^ = yhc/G explicitly, the Lorentzian action is (cf. [6| [16]. 
but note that here I am using A{t) for the physical scale factor and a{t) for the 
dimensionless rescaled scale factor, unlike the a{t) and r{t) used in [16] for those 
two respective quantities) 
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where b = a/3/ A is the radius of the throat of pure de Sitter with the same value 
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of the cosmological constant, A = A/(3m^) = l/(m6)^ is a dimensionless measure 
of the cosmological constant in units given by the mass of the inflaton, a = = 
mA = {uvf/^ = {UVf/^ and ^ = ^At^G/^cJ) = (1/4) In (t;/^) = (1/3) In (\//f/) 
are dimensionless forms of the scale factor and inflaton scalar field, u = e^""^*^ = 
^2g-2ip y _ g2Q!+2(/p _ ^2g+2¥3 gj.^ convenieut choice of null coordinates on the 



mmisuperspace (see, e.g.. 



U 



u 



3/4 



,(3/2)a-{3/2)<p 



i3/2e-2^ and V = v^l^ 



g(3/2)a+(3/2)ip _ ^3/2g+2(/j alternative choice of null coordinates, an overdot 

represents a derivative with respect to t, 5*0 = (37r)/(2G'm2) = (37r/2)(mpi/m)^, the 
DeWitt metric |20] on the minisuperspace is 
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rfs^ = ^oe^°(-c/a2^ V) = - 
the 'potential' on the minisuperspace is 

V = ^oe'"(</^' + A - e~2") = -So{uvy/^ 



--SodUdV, 
9 



(3) 



1 

16 



A + ^ln^^ 



u 



(4) 



alternative rescaled lapse functions are n = 2{uvY^^n = 2an 
nV = mNV, and the conformal minisuperspace metric is 



2m? AN and n 
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(5) 



The null coordinates u and v are chosen so that as one approaches the null boundaries 
of the minisuperspace, at m > 0, f = where a = with ip = — oo (except at 
u = V = 0, where ip can have any value), and at m = 0, f > where a = with 
p = +00 (again except at m = f = 0), the conformal minisuperspace metric Eq. (|5]) 
approaches {1 / 4) SQdudv , so that u and v are 2/5*0 times null coordinates that are 
the local analogues of orthonormal Minkowski coordinates near the boundaries. 

The Hamiltonian constraint equation and independent equation of motion can 
now be written as 



(j_dAy 

\Jla~dt ) 
l_d_ (]_dip 
Ndi \N~dt 



'I dpY 



,2_2 




+ m^p^ = 0, 



for general lapse function from the second form of the action above. 



p + 3-p + p = 0, 
a 



(6) 



(7) 
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from the third form of the action with n = 1, which will henceforth be assumed unless 
otherwise indicated (e.g., by including the lapse N or n exphcitly in a formula), 



• 2 -2 2 I \ 
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(p + 3d0 + = 0, 
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(8) 



for the fourth form of the action, and 
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for the minisuperspace null coordinates u = e^" and v = e^""'"^'^, and the alter- 
native null coordinates U = u^l^ = e(3/2)"-(3/2)^ ^nd V = v^l^ = e(3/2)°+(3/2)v_ xhis 
is with unit value for the rescaled lapse function, n = 1, but the last two equations 
appear simpler directly in terms of u and v if for just these equations we use the 
rescaled lapse n = 1 or n = l/(2a) = (1/2)6"" = {l/2){uv)-^/^ = {1/2){UV)-^/^, 
which gives 



uv 
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Although they are redundant equations, one may readily derive from Eqs. ([7]) 
and dHD that 

/d^ + l 



a{^^ - 2(p^ + A) 



20' 



and 



a 



-2a 



30^ 



11 



(12) 



when n = 1. Then when neither side of the constraint (first) equation part of Eqs. 
([HD vanishes (e.g., when V ^ 0), and when 7^ 0, one may define /' = df /dip = f /ip 
and reduce Eqs. ([H]) to the single second-order differential equation (cf. [6j) 



[a 



a 



l)((^a' + 3(/?2 + 3A-2e 



-2a\ 



ip^ + X 



-2a 



(13) 



Alternatively, when V ^ (or equivalently ^ ip'^), but when a ^ instead of 
7^ 0, one can write 



rfV {dip/daf - 1 



da'^ 



ip^ + X 



-2a 



3(^^ + 3A - 2e 



-2a 



dip 
da 



(14) 
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Yet another way to get the equations of motion is to note that the seventh 
(penuhimate) form of the action from Eq. ([2]) gives the trajectories of a particle 
of mass-squared V in the DeWitt minisuperspace metric ds'^, and the eighth and 
final form of the action gives timelike geodesies in the conformal minisuperspace 
metric ds"^ = Vds'^. When one goes to the gauge fi = 1, then [ds/dtY = — 1, 
so that along the classical timelike geodesies of ds^, the Lorentzian action is S* = 
— J dt = — J ^/—dP, minus the proper time along the timelike geodesic of ds^. 
However, one must note that the conformal metric ds^ = Vds^ is singular at = 
0, that is at <y9^ + A = e~^" = 1/a? = l/(mA)^, whereas there is no singularity 
in the DeWitt metric ds"^ or the spacetime metric along this hypersurface (curve) 
in the two-dimensional minisuperspace [a, (p) under consideration. The second- 
order differential equations ( IT3ll and ( IT4l) also break down at = and must 
be supplemented by the continuity of a and of (in a gauge in which n ^ is 
continuous there) across the V = hypersurface (curve). 

To get reasonable numbers for the dimensionless constants in these equations, I 
shall follow [H] and set m ^ 1.5 x IQ-^G-^/^ ^ 7.5 x 10-^{87tG)-^/^ [13 H], so the 
prefactor of the action becomes Sq = (37r)/(2Gm^) = (37r/2)(mpi/m)^ ^ 2 x 10^^, 
and the dimensionless measure of the cosmological constant is A = A/(3m^) = 
l/{mbY ~ 5 X 10~^^^. Thus A may be taken to be extremely tiny. 

The constrained Hamiltonian system for this ^ = -1-1 FLRWA-c/) model universe 
has an unconstrained 2 dimensional phase space with d = 2 that may be labeled 
by the coordinates and conjugate momenta Pi, which may be chosen to be any 
of the following sets: 

r. , Stt AdA 2n^A^d^^ 

a da dip 
[a, if, Pa = -i^o--j7^ = +'^0— -rrl, (16) 
n dt n dt 

{a, ip, Pa = -So — a, p^ = +So — 0}, (17) 
n n 

{u, V, pu = -—V, p^ = -—u\, (18) 
An An 

{t,, V, . J^v^ . J^U). (19) 
The Hamiltonian on this phase space is then 

3,^2 




AGm"^ 



-p^+pD + yne^"(-e-^" + A + v.^) 



6 



4 _ So_ 

Do 4 



A + -^In^ - I vuv 
lb u 



1 - A + - In 



V 



9 U 



(UV) 



2/3 



2 " 



n 



+ n (-a + Xa^ + a^v^^) 



-6one 



1 

4^ 



3a 



nj \n, 



+ A + 



-5*0 <J — n ^Mf — n 



A + ^ In^ - ) 



^So{-Vif7r-n(f/V)i/3 



u 



(20) 



The last four expressions are not in the canonical form as functions of the coor- 
dinates and momenta but are given to express the value of the Hamiltonian in terms 
of the coordinates and their time derivatives. Because the Hamiltonian constraint 
equation, obtained by varying the action S of Eq. ([2]) with respect to the lapse 
function n or n, is if = 0, these last four expressions for H can be easily seen to 
lead to the first equations in each of Eqs. (I7j), (jH]), (E]), and ( fTOl) when one chooses 
the rescaled lapse function n to be 1 for Eqs. ([7]), ([8]), and ([9]) and chooses n = 1 in 
Eq. (fTOj) . One can also write the Hamiltonian constraint equation H = directly n 
terms of the canonical coordinates and momenta as 



22 2 o2 6 f \ , 2 
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(21) 
(22) 
(23) 

(24) 



3 The canonical measure for /c = +1 FLRWA- 



The natural symplectic structure for the k = +1 FLRWA-0 constrained Hamiltonian 
system is the closed nondegenerate differential 2-form 

= ^2 = dPi A dQ^ 

= dpA AdA + dp^ Ad(f)= -Som^AdA AdA + 27i'^md{A^^) A dcj) 

= dpa A da + dp^ A d(p = Sq (^—add A da + d{a^ip) A dipj 

= dpa Ada + dp^ Adip = Soe^^i^—da A da + dip A dip + 3ipda A dip) 



dpu A du + dpy A dv = —-^{uv 



a/4 



dv A du + dii A dv + - { ] du A dv 

A\v uJ 
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= dpu A dU + dpv A dV = --So [dV A dU + dU A dV) , (25) 

where for the expressions in terms of the time derivatives, I have used the default 
option n = 1. 

When this is pulled back to the H = constraint hypersurface of dimension 
2(i — 1 = 3 in the unconstrained phase space Td = r2 of dimension 2d = 4 and fur- 
ther pulled back to an initial-data surface r^^.i = Fi of dimension 2d — 2 = 2 that is 
transverse to the Hamiltonian flow in the 3-dimensional constraint hypersurface, it 
gives the symplectic structure differential form u = Ud-i = oJi on that 2-dimensional 
initial data surface. Since li — 1 = 1, it is the first power of this symplectic structure 
form that gives the canonical Liouville-Henneaux-Gibbons-Hawking-Stewart mea- 
sure or volume (area) element ^n-i = f^i = w on the initial data surface Fi [2]. 
That is, if a bunch of orbits B intersects an initial-data surface S in the region S, 
the canonical measure of that bunch of orbits is /^(-B) = Js^- As Gibbons, Hawk- 
ing, and Stewart show [2], this measure is preserved as one follows the bunch B to 
where it intersects a different initial-data surface S' in the region S', giving the same 
measure yu(-B) = Js'<^- 

Here we are restricting to nonsingular cosmologies, Friedmann-Lemaitre- Robertson- 
Walker universes that have neither a big bang nor a big crunch, so the scale factor 
A or a never goes to zero. Except for a discrete set of zero canonical measure [6], all 
of these solutions will contract from infinite a at infinite past time and re-expand to 
infinite a at infinite future time. Therefore, they will each have a global minimum 
for a = e", which I shall label = exp {am), where da/dt = d = e"d = and 
hence Pa = Pa = Pa = 0. Let ipm, 0m, «m, and be the values of ip, ip, a, and 
at this global minimum for a and hence also for a. 

The Hamiltonian constraint H = 0, given by the first Eq. ([7]) with both sides 
equal to zero, implies that (using the default setting of the rescaled lapse function 
as n = 1) 

a^ = (<^^ + 0^ + A)-^/^ (26) 

or 

am = --\n{^l + ipl + X), (27) 

Therefore, initial data are given by values of and with = dmi'Pm^'Pm) 
given by Eq. fl26|) . and then one has a = automatically at this point in the 
constrained phase space. 

One can further readily see (cf. [I2l [T3])that at an extremum for a, where d = 0, 
that one has 

- = a = 3(^2 ^ 3A - 2e-2" = - 20^ + A. (28) 
a 

Therefore, for an extremum to be at least a local minimum, one needs < 
(2/3)^/^((y9^ + A)"^/^ or ipl^ < {ifl^ + A)/2, though there are additional conditions 
for such a local minimum to be a global minimum. 
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Since Eq. (!26|) gives a unique value for am for each set of real values for (pm 
and ifrn, it naively appears that there are no constraints on and (fm- However, 
different sets of and ipm can lead to the same solution, since a solution may 
have more than one point along its trajectory (more than one time) where a = 0, 
only one of which may be a global minimum in the generic case in which there are 
not more than one time with the same global minimum value of a{t). Therefore, 
counting all possibilities for ipm and ipm overcounts the trajectories that have d = 
somewhere along them. Furthermore, singular trajectories, with a going to zero in 
the past or future, may also have points where d = and thus be counted if one 
counts all possible pairs of (fm and (fm- In the next Section we shall look at the 
restrictions on and ipm in order that this pair correspond to a global minimum 
of a rather than some other extremum like a local minimum that is not a global 
minimum, or either a local or a global maximum. However, first we shall show that 
the total canonical measure of all solutions with any local nonzero extremum for a 
(at a = am or a = am where d = and a = 0) has finite measure. 

It is convenient to define two new variables (3 and 6 so that (with rescaled lapse 
function n = 1 as usual) 



If = e-^sinO. (29) 
The Hamiltonian constraint equation Eq. ([8]) then becomes 

a^ = X + - (30) 
One can also easily calculate that the time derivatives of /3 and 6 are 

/3 = 3asin^9 = ^a{l-cos2e), (31) 

^ = -1 - 3dsin^cos^ = -1 - ^dsin2^. (32) 

When averaged over one period of ^ in a regime in which the scalar field os- 
cillates rapidly relative to the expansion (so that the time-average of the scalar 
field stress-energy tensor is approximately that of pressureless dust), /3 changes by 
approximately 3/2 as much as a, so it is convenient to define a total rationalized 
dimensionless 'mass' (twice the energy density multiplied by the volume and divided 
by 5*0 and by the scalar field mass m) that is nearly constant in the dustlike regime 
il6j: 

M = e'''-'^ = a%ip' + ip'), (33) 

obeying 

M = 3Mdcos2^ (34) 
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or 

^ = 3cos2^. (35) 
da 

The symplectic structure 2-form Un = ^2 given by Eq. fl25|) is written in terms of 
the four independent 1-forms of the unconstrained phase space r„ = r2 of dimension 
2n = 4. When one imposes the Hamiltonian constraint H = 0, one of the four 1- 
forms appearing in U2 can be written in terms of the other three. Thus one can write 
the symplectic structure 2-form in terms of three 1-forms that are independent on the 
constraint hypersurface H = 0. Choosing these three to be various combinations 
of the differentials of a, a, ip, ip, /3, 6, and M, one can write the 2-form on the 
constraint hypersurface as 

cu = Soe^°' ( -^da A dip + ^da A dip - dip A d(p] 
\ da da J 

= Sqc^" {—ipda A dip + Lpda Adip — adip A dip) 

= Soe^""-^^ ( ^da A d^ - ^da A de + d(3 A de] 
\da da J 

= Soe^""-^^ (eda A d(3 - pda A dO + ad(3 A dO) 

= -Sol-^daAdM + ^daAde-dMAde]. (36) 
2 \ da da J 

Here the derivatives with respect to time (with = 1) or to a that are the coefficients 
of the basis 2-forms inside the parentheses are derivatives along the trajectories 
forming the cosmological spacetime solutions, unlike the basis 1-forms that make up 
the basis 2-forms that are differentials transverse to the trajectories. 

On an initial data surface that is an extremum of the scale size or of the logarithm 
of the rescaled scale size, a = am, where a = 0, one has /3 = /3m, = 9m, and 
M = Mm = e^"™"^^", and from the two coordinates /3m and 9m on this initial 
data surface that I shall call Se, one gets am = —(1/2) In (A — e^^'^™) or /3m = 
-(l/2)ln(e-2°'" - A), then giving Mm = e'''"(l + Ae^^'")-^/^ = e"'"(l - Ae^"™). 
Note that for an extremum we must have am < —(1/2) In A or am < 1/"\/A or 
Aa^ < 1, but /3m can be an arbitrary real number (though one no longer has the 
full range of all real numbers for /3 if the extremum is required to be a global 
minimum for a or a). One can alternatively label the initial data surface Se of 
extrema (all points in the constrained hypersurface H = where also d = 0) by the 
two coordinates am and 9m, both of which have finite ranges, < < 1/"\/A and 
< 9m < 2-iT. Then on that initial data surface Se one has = In am, 

/3 = /3m = -(l/2)ln(l/a^- A), (p = (pm = e'^^ cos 9m, ^ = (pm = e-^"'sm9m, 
am=[l- 3(1 - XalP)sm'^9m]/am, and Mm = ^^(1 - Aa^). 

The pullback of U2 to an initial data surface Se that is at an extremum of a 
and of a, where d = and hence where am = ami^m, P'm) = (A + 0^ + v^^)"^^^ = 
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(A 



< 1/ A i 



IS 



Soe"^--^^-a^dl3^ A de^ 
- 3 f 1 - Aail sin^ e„ 



Sq 



3 1 -x^ sin^ a 



dttm A dOn 
dx A d9m, 



(37) 



where /io = So/y/X ^ 3 x 10^"^ and x = \/Xam, which has the range < a; < 1. 

Because both (or x) and 6m have finite ranges, and because the integrand is 
bounded above within this range, the measure for the set of solutions with a nonzero 
extremum for a is finite. For the case with zero cosmo logical constant. Hawking and 
I showed [3J that the solutions with an extremum within a finite range of a have 
finite measure, but for that model there is no upper bound on a at an extremum, 
and almost all solutions have a maximum for a, so the total measure for solutions 
with maxima is infinite. But in the present case, the positive cosmological constant 
imposes an upper limit on the value of a at an extremum. 

Therefore, we see that the set of solutions with a nonzero extremum for a has a 
finite canonical measure, out of the infinite measure for all solutions for a k = +1 
FLRW cosmology with a minimally coupled massive scalar field and a positive cos- 
mological constant. The finite measure of solutions includes both totally nonsingular 
solutions, which have a nonzero global minimum for a, and also solutions with both 
a big bang and a big crunch, which have a finite global maximum for a. It also 
includes solutions that start at a big bang and eventually expand forever, and their 
time reverses that contract from a = oo to a big crunch, so long as they have a local 
extremum for a, where a or the Hubble variable a is zero. The only set of solutions 
that have infinite measure are those that expand monotonically from a big bang at 
a = to infinite size at a = oo, or the time reverses that contract monotonically 
from a = cxD to a big crunch at a = 0. 



4 Canonical measure for nonsingular cosmologies 

We have seen that the total Liouville-Henneaux-Gibbons-Hawking-Stewart canoni- 
cal measure for nonsingular Friedmann-Lemaitre- Robertson- Walker cosmologies with 
a minimally coupled massive scalar field is finite. (For there to be a nonzero measure 
for such nonsingular cosmological solutions of the Einstein-scalar field equations, we 
need a closed cosmology with = -|-1 to allow a to have a minimum value, and 
we need a positive cosmological constant to allow a to go to infinity asymptotically 
in both directions of time.) Now let us calculate the measure for the nonsingular 
solutions. 

The canonical Liouville-Henneaux-Gibbons-Hawking-Stewart measure or volume 
(area) element u on the initial data surface Fi with a = 1, given by Eq. f l37p . has 
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a sign of the coefficient of dx A dOm that is proportional to a = cim = [1 — 3(1 — 
Aa^) sin^ Om\/am = V^[l — 3(1 — x^) sin^ the acceleration of the scale factor 

a at its extremum. When this is positive, the extremum is a local minimum for the 
scale factor; when dm < 0, the extremum is a local maximum for a. If we integrate 
oj over the range giving positive dm, we get the finite measure /xi = {An / \/T{)iiq f» 
2.4184yUo. If we reverse the sign of u and integrate it over the range giving negative 
dm-, we get the same finite measure, /i2 = (47r/-\/27)/io- Therefore, the total measure 
for solutions with extrema for a is not greater than (13=^1+^2 = (87r/-\/27)/^o, 
finite. 

However, solutions may have more than one extremum for a, and /is counts all 
such solutions with a multiplicity given by the number of extrema that they have. 
Therefore, let us calculate what the measure is for nonsingular solutions by just 
taking the measure at the nonzero global minimum for a, avoiding the overcounting 
of a multiplicity of minima. For this calculation we shall assume that A ^ 1, as it 
indeed is in our part of the universe where A = A/(3m^) ^ 5 x 10~^^^, and hence 
drop correction terms proportional to positive powers of A. 

Most of the measure for the nonsingular solutions will come from values of a not 
too much less than the maximum value for an extremum, which is at a = 1/ a/A or 
X = y/Xttm = 1- Therefore, we can assume that x is not enormously smaller than 
unity for almost all of the measure. For a nonsingular solution that has a global 
minimum at a = Ina = am = Inx — (1/2) In A ^ 1 (with —(1/2) In A ~ 127 ^ 1), 
the fact that Eq. ( l3Ti) implies that /3 cannot decrease as a increases implies that the 
Hamiltonian constraint equation Eq. ( l30l) gives 

a^ = X + _ < A + e-'^™ - e"'" = e-^"™ - e"'" < 4- (38) 

x^ 

For X > ^ 7 X 10"^*^, we thus get < 1. 

As a result, Eq. (!32l) implies that ^ — 1 to high accuracy, and Eq. implies 
that M = e^"~^^ = 0^(9?^ + ip^) stays very nearly constant along most of the 
nonsingular trajectories. Eq. fl3T|) implies that P is also very small, so over a number 
of oscillations of the scalar field that is not too large (a change in the phase 6 that 
is not too many times 27r), neither a nor /3 change much, though after a very long 
time and a huge number of oscillations of the scalar field (enormous change in 6), 
both a and (3 grow indefinitely, while M and %p = 9 + t stay nearly constant and 
indeed both approach precise constants in the infinite future, M^o and ipoo- 

(To define ipao unambiguously, set t = at the global minimum for a and require 
< 6*™. < 27r there. One can make this definition not only for nonsingular solutions 
but also for big bang solutions that start at global minimum for a that is a = 1, 
where one can set t = 0, and then evolve to infinite a where ipoo can be evaluated. 
To circumvent the jumps in 9m at the minimum that would occur with a sequence 
of solutions with 9m approaching 27r and then jumping back to 0, instead of defining 
the two real constants M^o and ipoo-, it would be better to define the one complex 
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constant Z = a/MooC*'^"", which is invariant under shifting 6^ and hence ip and ipoo 
by an integer multiple of 27r. Any solution that evolves to a = oo will have a definite 
value for Z that may be obtained by analytic integration of the equations of motion 
from any initial point in the constrained phase space, except for the nonsingular 
solutions that have two equal global minima for a and therefore the ambiguity of 
which one to use for setting the zero of t, and the solutions that are the limits of 
sequences of solutions with bounces of arbitrarily large values of the scalar field [6]. 
Both of these types of particular solutions will have zero measure, so all but a set 
of measure zero of the solutions that evolve to a = oo will be integrable, having 
two real constants of motion that may be given by one complex constant Z, that 
are analytic functions over all but a set of hypersurfaces of measure zero in the 
constrained phase space. The same will be true for k = and k = —1 FLRW-scalar 
models with a nonnegative cosmological constant to allow solutions to expand to 
a = oo, though since they cannot have extrema of a that in the k = +1 case can 
lead to hypersurfaces of the constrained phase space where the constants of motion 
are not analytic, it appears that the k = and k = —1 FLRW-scalar models will 
be totally integrable over the entire constrained phase space. In fact, since these 
models give a expanding monotonically from a = to a = oo, or the time reverse, 
one can not only define constants of motion by the asymptotic behavior of M and 
ip at a = oo that gives rise to the complex constant Z, but also by the asymptotic 
behavior at a = 0, such as the value ofv—u where either one of these null coordinates 
u and V goes to zero, and the value of the slope dv/du there.) 

During the oscillations of the scalar field while a and /3 stay near their values at 
the extremum, one can write 



d'^a d'^a 



de^ 

Integrating this gives 



where 



and then finally 



dt^ 
da 

de 



3x' 



2x2 
2^ 



B 



3(1 -x2) 
2(3x2 - 1)^ 



l + 3(l-x^)(l-cos2^) 



[9 - 9o + B sm29)] , 



9o = 9m + Bsm29m, 



a 



A(3x2 - 1) 
4x2 



Bcos29 - (9n-9r. 



B cos 29n 



(39) 

(40) 
(41) 

(42) 



One can then see that for = to be not only a local minimum (which 
requires merely 3x2 _ ^ _^ _ ^2^^-^ _ cos26'm) > 0) but also a global minimum, 
one needs that 3x2 _ i y q ^qj. nonnegative coefficient of the quadratic term in 9. 
Furthermore, by sketching the behavior of a{9), one can see that for — vr < 29q < it, 
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one needs — 7r/2 < 20m < 7r/2; for vr < 26q < Svr, one needs 37r/2 < 26m < 57r/2; for 
Stt < 26q < Stt, one needs 77r/2 < 26*^ < 97r/2; etc. Thus we cannot have a global 
minimum with tx/2 < 29m < 37r/2, 57r/2 < 2^^ < 77r/2, etc. 

If we choose 6m to lie in the range < < 27r, then there are four allowed 
ranges of 6m of width tc/A (covering half the full circle for 6m', the other half does not 
give extrema that are global minima) that give equal contributions to the measure. 
Let us focus on the first, which is that part of < 6'm < 7r/4 that gives 6o = 
6m + Bsm26m < 7r/2. Using Eq. fHTl) to express B in terms of x\/Xam allows one 
to convert this to a restriction on x for < < 7r/4: 



Xm{6m) = \ ^ 77^ 1 o ■ o/j < < 1. 43 

V 3vr — 66m + 3 sm 26m 

If we now integrate the canonical Liouville-Henneaux-Gibbons-Hawking-Stewart 
measure or volume (area) element u in Eq. ( 137|) over the initial data surface, say 
Sm, that has a not only an extremum but also a global minimum (which is 4 times 
the integral of u over the one region above, in order to include all possibilities for 
< 6m < 277 which give a global minimum), we get (using y = 26m) 



u} = fiol2-J^ dy {3 cos y - l)xm + {I - cos y)Xm\ . (44) 



Doing the integral numerically with Maple 12 gave 7 ~ 0.86334. Putting in 
the numbers given above for 5*0 = (37r)/(2G'm^) = (37r/2)(mpi/m)^ 2 x 10"*^^, 
A = A/(3m2) = l/(m6)2 ^ 5 x 10"^", and /io = Sq/VX ^ 3 x 10*^^ gives the 
measure for the nonsingular k = +1 FLRW cosmologies with the observed value of 
the cosmological constant and a scalar field mass of m ~ 1.5 x 10~^mpi as 

S 3 /3 

l^m = 7/^0 = 7^ = 7 ,/ ~ 0.86334/io ^ 3 x 10^^ (45) 

To convert this to a number of quanta, say Nm, one divides the phase space measure 
hj h = 2Txh = 27r in our units with h = c = 1 to get 

iV. - - - 7-% - 7-^ - 0.13740/.0 - 4 x 10^. (46) 
27r 27r 27rVA AGmy h. 

One can compare this with the maximum number of scalar field quanta, say 
Nm, that one can have for a nonsingular = +1 FLRW cosmology with a positive 
cosmological constant if the scalar field acted as pressureless dust, which is how it 
does act at late times when one averages over an integer number of oscillations of the 
scalar field. Then M would stay constant. For a universe with a minimum value of a 
that is ttm, one gets M = 0^(1 — Aa,^), which has a maximum value (when 1— 3Aa,^ = 
l-3x^ = 0) of MAf = 2/yi7A. The physical energy density (l/2)[(m0)2 + (0/A^)2] = 
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[3/(87rG)]m^((^^ + LfP') (with N = n/m = 1/m) multiplied by the physical 3- volume 
27r^A^ = then gives a physical 'mass' M. — {So/2)mM, so the maximum 

number of scalar dust quanta of mass m in a nonsingular A; = +1 FLRW cosmology 
with a positive cosmological constant is 

Nm = 7:SoMm = = ~ 0.19245/Xo ^ 1.400607Ar„ ^6x 10^^ (47) 

2 v27A v27 

Therefore, the number of quanta corresponding to the actual phase space measure 
over the nonsingular k = +1 FLRW cosmologies is \fTt^l (27r) 0.713976 times the 
maximum of that for a dust model with the same particle mass. In the actual scalar 
field model, the scalar field undergoes coherent oscillations in which the phase has 
gravitational consequences, so it cannot be accurately modeled by assuming that 
the scalar field is in a precise number eigenstatc with a totally undetermined phase, 
which would give zero pressure for a homogeneous field such as is being assumed 
here. 



5 Canonical measure for inflationary cosmologies 

Nearly all of the finite total Liouville-Henneaux-Gibbons-Hawking-Stewart canonical 
measure for nonsingular Priedmann-Lemaitre- Robertson- Walker cosmologies with a 
minimally coupled massive scalar field and a positive cosmological constant occurs 
for solutions that are not large deviations from empty de Sitter spacetime. For 
example, if one defines the rationalized dimensionless energy density to be 

+ (48) 

which is (87rG)/(3m^) times the physical energy density p = (l/2)[(m0)^-|-(0/A^)^] = 
[3/ (87rG)]m^(99^ + 0^) = [3/ (87rG)]m^p of the scalar field, and sets p = Pm = o-m^ — A 
at the global minimum for the scale factor at a = a^, then under the approximation 
that M is constant, one has l/\/3A < < 1/\/A and p < Pm ^ '^^ everywhere in 
the spacetime, so the physical energy density of the scalar field is never more than 
twice the physical energy density pA — A/(87rG) = [3/(87rG)]m^A corresponding to 
the cosmological constant, that is. p < 2p\. 

If Pm = — A > 2A ^ I at an cxtremum of a, then for pm, <^ Ij ^ 1 for the 
entire solution, so M stays nearly constant at a value less than its maximum value 
for nonsingular dust solutions [M = Mm = 2/V27A), and the resulting solutions 
collapse to a = rather than expanding to infinity. To obtain solutions with p^^ > 
2pA that expand to infinity in both directions of time rather than collapsing to zero 
size, one need to have a period of inflation in which M grows larger to become larger 
than Mm = 2/V27A, as one can see from the following argument: 
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The Hamiltonian constraint equation in Eq. ([7]) can be written in terms of a and 
M as 

d2 = /(a) = Aa^ - 1 + — . (49) 



a 

> 



After the end of a possible period of inflation (which requires p ~ 1 ^ A), during 
which M can grow exponentially, M will become nearly constant as the scalar field 
starts oscillating with a period much less than the inverse of the Hubble expansion 
rate. Then the universe will expand forever if /(a) stays positive for all larger a. 
For constant M, the extremum of /(a) is at a = [M/(2A)]^/^ = (M/MmY^^/V^- 
The value at the extremum is /(a) = (e.rSM^A)^/^ - 1 = {M/MmY'^ - 1, so if 
inflation ends before a reaches the extremum, one needs M > Mm in order for 
d to stay positive as a passes through the extremum of /(a). If inflation gives 
M < Mm, it will end far before this value of a is reached, so one needs inflation to 
give M > Mm = 2/y/27X > 1 if one starts at an extremum with > 2A and hence 
with < 1/ "\/3A, where there are no noninflationary solutions (M nearly constant) 
that expand to infinity in both directions of time and hence are nonsingular. 

For a symmetric bounce {ipm = or 9m = 0) at ipm = v^b ^ 1, I have calculated 
numerically [16] that the number of e-folds of inflation (not be be confused with 
the previous use of for the lapse function) is 

N{^t) + \ In - 1.0653 - - ^, (50) 

and that the asymptotic value of M is 

To give M > Mm, this requires A^ > Nm ~ 44.27. 

For a = A+pm(l— 3 sin^ 9m) ~ pm(l— 3 sin^ 9m) > 0, we need | sin6'm| < l/y/S. If 
this is is true for pm^ I and 9m is not too close to the boundary, then the equations 
of motion in this inflationary regime will generally cause (p and hence 9 to decay to 
near zero, and then one will get roughly A^ ~ {3/2)pm e-folds of inflation, though 
the coefficient (3/2) will become some 6'm-dependent number that is (3/2) only at 
sin^m = 0. However, for a crude estimate of the measure for different amounts of 
inflation, let us ignore this effect. Then Eq. f l37|) with = (A + Pm)^^^^ ~ Pm^^^ 
gives the measure of varying numbers A^ of e-folds of inflation as 



If 

--5'o j{X + Pm)~^^'^ A + Pm (1 - 3sin^ 6'rrt) dpmAdifr 

Soil-3smHm)diPm'^^)Adipm 
(y6-cos-^3-^/^)5o/d(p„V2) 
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(3 - 73/2 008-^ 3-^/^)So J d{N-'/^). 

SojdiN-'/'). (52) 



The second expression on the right hand side (after / u) uses the very good approx- 
imation that \ <^ for the inflationary values of Pm that are at least of the order 
of unity in our units that set effectively set m = 1. The third evaluates the integral 
over 9jn- The fourth uses the approximation N (3/2)p^ that is actually only true 
for sin = 0, so the fifth drops the uncertain numerical coefficient. 

Thus we see that the measure for at least e-folds of inflation is proportional 
to 1/ for large N . If we take the fraction of the total measure for nonsingular 
solutions, which was "^fp^ = ^So/VX a; 3 x 10®^, the fraction of the total measure 
for at least e-folds of inflation is 

^ _ yxT^ _ 10-^^ 10-56 ^^^^ 



That is, the fraction of the measure for all nonsingular k = +1 FLRW cosmologies 
that have inflation (requiring N > Nm ~ 44.27 for the present toy model with just 
the observed value of the cosmological constant and a massive scalar field with 
m ~ 1.5 X IQ-^mpi; it would be higher if the scalar field could decay into radiation 
at the end of inflation) is about 10^''"^'^. However, the fraction goes down with the 
minimum number of e-folds required, N, only by an inverse square root of A^, and 
not as e"^^, so there is no conflict with not observing the universe to have such a 
minimal amount of inflation that spatial curvature is observable. 



6 Conclusions 

Although the total canonical Liouville-Henneaux-Gibbons-Hawking-Stewart mea- 
sure is inflnite for Fricdmann-Lemaitre-Robertson- Walker classical universes with a 
minimally coupled massive scalar fleld and a positive cosmological constant, it is 
flnite for completely nonsingular solutions (which have positive scale factor every- 
where) . Nearly all of the solutions have the energy density never more than twice the 
effective energy density of the cosmological constant, but the tiny fraction, ~ lO^^^, 
of the measure in which the energy density ever exceeds this tiny amount has at 
least ~ 44 e-folds of inflation and gives a measure that decreases only very slowly 
(as an inverse square root) with the minimal number of e-folds required. 



17 



7 Acknowledgments 



I am grateful for the hospitality of Princeton University, where I was motivated 
to work on this problem by a talk by Neil Turok on his latest approach to the 
measure issue [5]. This research was supported in part by the Natural Sciences and 
Engineering Research Council of Canada. 



References 

M. Henneaux, "The Gibbs Entropy Production in General Relativity," Nuovo 
Cim. Lett. 38, 609-614 (1983). 

G. W. Gibbons, S. W. Hawking, and J. M. Stewart, "A Natural Measure on 
the Set of All Universes," Nucl. Phys. B 281, 736-751 (1987). 

D. N. Page and S. W. Hawking, "How Probable Is Inflation?" Nucl. Phys. B 
298, 789-809 (1988). 

G. W. Gibbons and N. Turok, "Measure Problem in Cosmology," Phys. Rev. 
D 77, 063516 (2008). 

N. Turok, in preparation (2011). 

D. N. Page, "A Fractal Set of Perpetually Bouncing Universes?" Class. Quant. 
Grav. 1, 417-427 (1984). 

L. Parker and S. A. Fulling, "Quantized Matter Fields and the Avoidance of 
Singularities in General Relativity," Phys. Rev. D 7, 2357-2374 (1973). 

A. A. Starobinsky, "On a Nonsingular Isotropic Cosmological Model," Sov. 
Astron. Lett. 4, 82 (1978). 

S. W. Hawking, "Quantum Cosmology," in Relativity, Groups and Topology 
II, Les Houches, 1983, edited by B. S. DeWitt and R. Stora (North- Holland, 
Amsterdam, 1984), pp. 333-379. 

V. A. Belinsky, L. P. Grishchuk, I. M. Khalatnikov, and Ya. B. Zeldovich, 
"Inflationary Stages In Cosmological Models With a Scalar Field," Phys. Lett. 
B 155, 232-236 (1985). 

V. A. Belinsky and I. M. Khalatnikov, "On the Degree of Generality of In- 
flationary Solutions in Cosmological Models with a Scalar Field," Sov. Phys. 
JETP 66, 441 (1987). 

A. Yu. Kamenshchik, I. M. Khalatnikov, and A. V. Toporensky, "Simplest 
Cosmological Model with the Scalar Field," Int. J. Mod. Phys. D6, 673-692 
(1997). 



18 



A. Yu. Kamenshchik, I. M. Khalatnikov, and A. V. Toporensky, "Simplest Cos- 
mological Model with the Scalar Field II. Influence of Cosmological Constant" 
Int. J. Mod. Phys. D7, 129-138 (1998). 

N. J. Cornish and E. P. S. Shellard, "Chaos in Quantum Cosmology," Phys. 
Rev. Lett. 81, 3571-3574 (1998), |gr^qc79708046, 

A. Yu. Kamenshchik, I. M. Khalatnikov, S. V. Savchenko, and A. V. Toporen- 
sky, "Topological Entropy for Some Isotropic Cosmological Models," Phys. Rev. 
D 59, 123516 (1999) jir^qc79809048| 

D. N. Page, "Symmetric-Bounce Quantum State of the Universe," JCAP 0909, 
026 (2009) larXiv: 0907. 1893] . 

A. Linde, Particle Physics and Inflationary Cosmology (Harwood Academic 
Publishers, Chur, Switzerland, 1990). 

A. R. Liddle and D. H. Lyth, Cosmological Inflation and Large-Scale Structure 
(Cambridge University Press, Cambridge, 2000). 

D. N. Page, "Minisuperspaces with Conformally and Minimally Coupled Scalar 
Fields," J. Math. Phys. 32, 3427-3438 (1991). 

B. S. DeWitt, "Quantum Theory of Gravity. 1. The Canonical Theory," Phys. 
Rev. 160, 1113-1148 (1967). 



19 



